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model for a finite region of the parameter space surrounding the critical 



where RKKY interaction is comparable to the Kondo temperature. After 
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, several exact simplifications involving reduction to one dimension and abelian 
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bosonization, the critical point is explicitly identified, making clear its physi 
cal origin. By using controlled low energy projection, an effective Hamiltonian 
is mapped out for the finite region in the phase diagram around the critical 
point. The completeness of the effective Hamiltonian is rigorously proved 
from general symmetry considerations. The effective Hamiltonian is solved 
exactly not only at the critical point but also for the surrounding Fermi-liquid 
phase. Analytic crossover functions from the critical to Fermi-liquid behavior 
are derived for the specific heat and staggered susceptibility. It is shown that 
applying a uniform magnetic field has negligible effect on the critical behav- 
ior. A detailed comparison is made with the numerical renormalization group 
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and conformal field theory results. The excellent agreement is exploited to 
argue for the universality of the critical point, which in turn implies universal 
behavior everywhere inside our solution region. 
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I. INTRODUCTION 



For a vast number of materials with strong electron correlation, the low energy excitations 
involve both itinerant electrons and well localized magnetic moments residing periodically on 
the lattice sites. This is the case of heavy fermion compounds JT] , and to certain extent it is 
also the case of high temperature superconducting cuprates . In such systems, two effects 
have crucial influence on the low energy properties and they compete with each other. They 
are the Kondo effect || and RKKY interaction M , which represent two different tendencies 
of the system to quench the local moments with conduction electrons or by themselves. The 
simplest model capturing both effects is the two-impurity Kondo model ||-(7|]. It is also 
believed that any possible new physics that may occur in lattice due to the competition 
between the two effects should be contained in the two-impurity problem ||. 

A simple way to see how the competition arises in the two-impurity Kondo model is to 
look at the problem from the scaling point of view || . The complexity of the problem will be 
fully realized this way and the concrete task we are facing will be defined. For all practical 
purposes, the bare Kondo coupling constant Jo and the RKKY interaction Kq are much 
smaller than the Fermi energy ep. Thus, if we form two dimensionless coupling constants 
with the help of the conduction electron density of states p_f(~ ^/ € f), PfJq and PfK , they 
are always in the weak coupling regime. However, as we start to eliminate high energy con- 
duction electron states near the top and bottom of the conduction band, both dimensionless 
coupling constants grow under renormalization and they mutually renormalize each other. 
Simple dimensional counting shows that pp Jo has dimension one and is marginally relevant, 
while PfK has dimension zero and is relevant. For these two relevant interactions we can 
define two energy scales, the Kondo temperature Tk and RKKY temperature Tricky, such 
that they correspond to the values of the decreasing conduction bandwidth at which the 
renormalized dimensionless coupling constants ppJ and pfK reach the unity respectively. 
In either Tk 3> Trkky or Tk <C Trkky limit, the problem is simple because we can perturb 
one of the two interactions. The most difficult situation corresponds to T K ~ T RK ky- This 
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is also the situation of most practical interest. 

Due to the broad interest in the competition between the Kondo effect and RKKY 
interaction, extensive investigations have been carried out in the last decade [|5yiO|-|T8|1 . As 
a result, a convincing phase diagram has emerged [pTT] , |T^| -p^1 , if not yet universally accepted 
without reservation This phase diagram is shown in Figure |TJ. The model exhibits 

Fermi-liquid behavior everywhere except at a special point on the particle-hole symmetric 
axis. At this point, the ratio between the fully renormalized effective RKKY interaction 
and the Kondo temperature is numerically estimated to be 2.2 [|ll]. The effective RKKY 



interaction is actually the RKKY temperature whose meaning we have explained in the 
above. But for some reason it has been confusingly called RKKY interaction. Although the 
precise numerical value of the ratio may depend on the individual's convention of defining 
the coupling constants, the important message is that this critical point is located at Tk ~ 
Trkky- 

Although several asymptotically exact results have been available about the critical point 
in certain limits [TT|,|I8|] , the exact physical origin of this critical point had not been unveiled 
until recently |2t| . By an explicit identification of two local impurity spin states whose level 
crossing being the origin of the critical point, we have rigorously shown how the constraints 
set by the discrete symmetries of the model ensures the occurrence of a non-Fermi-liquid 
critical point. We have also presented an effective Hamiltonian for the finite solution region 
as marked in Figure [I], and listed the low temperature properties of the critical point. In 
this paper, we present a detailed derivation of the effective Hamiltonian, and for the first 
time a full analytic solution for the whole solution region of Figure [|. 

Having admitted the existence of the critical point, we can already present a framework 
for the solution inside that solution region in Figure ^ by only invoking general scaling 
ideas. It is then the task of section [TV] to fill in concrete results. Since Kondo effect always 
takes place in our solution region, the basic energy scale must be the Kondo temperature 
Tk, which is much smaller than the Fermi energy. On the scale of Tk, the system has 
already lost its memory of microscopic details existed on the energy scale of the Fermi 



energy. The above mentioned Trkky is of the same order as Tk, and therefore does not 
constitute a new energy scale by itself. However, the competition induces the second energy 
scale T c . Inside our solution region, T c <C Tk, and T c vanishes at the critical point. For 
any physical quantity, its dependence on the bare parameters of the Hamiltonian should be 
absorbed into these two energy scales. For instance, we can write the specific heat in the form 
C(T) = f(T /Tk,T c /Tk), where f(x,y) is some universal two- variable function. The role of 
Tk is simply to set an energy unit for the problem. The second energy scale T c determines 
the crossover from the non-Fermi-liquid behavior governed by the critical point at T ^> T c 
to the Fermi-liquid behavior governed by the stable Fermi-liquid fixed point at T T c . 
If we recall how the appearance of the Kondo temperature T K (<^. ep) in the one-impurity 
Kondo problem leads to the drastic enhancement of various physical quantities including 
specific heat, we can expect additional enhancement from the appearance of T C (^C Tk)- 
The translation of this effect to the lattice problem will be a new mechanism for the heavy 
electron mass. 

Moving far away from the critical point in the phase diagram, the accuracy of our solution 
deteriorates. However, the low energy exponents for all physical quantities should not change 
since after all the system is still governed by a Fermi-liquid fixed point, as in the solution 
region of Figure [l] near the critical point. What need to be improved are the constant 
prefactors. Usually, a physical quantity calculated for a Fermi-liquid fixed point is a sum 
of several contributions with the same exponent. Among them, only a few are associated 
with the energy scale T c , while the others are associated with Tk- What we calculate in 
this paper are those contributions associated with the energy scale T c . This is enough near 
the critical point because they are enhanced inside our solution region of Figure [l|. As one 
moves away from the critical point, the other contributions associated with the energy scale 
Tk become increasingly significant. Certainly, one can always fit the prefactors for every 
physical quantity if well established numerical results or experimental data are available. 
But it is the merit of the theory to establish the relations between these prefactors in the 
same spirit of Nozieres's Fermi-liquid theory of Kondo effect [pi ]. A complete solution of 



the two-impurity Kondo model for the whole parameter space is beyond the scope of this 
paper. 

The layout of the paper is the following. In section [H], we present the preliminary trans- 
formations on the two-impurity Kondo model including the reduction to one dimension and 
bosonization. In section [TTT| , we identify the critical point and derive the effective Hamil- 
tonian. A rigorous proof of the completeness of the effective Hamiltonian is included. In 
section |TV| , we solve the effective Hamiltonian and calculate the low energy thermodynamics. 
In section [V], we compare our results with those derived from the numerical renormalization 
group and conformal field theory approaches. The universality of the critical point is strongly 
argued for. We conclude our paper in section [Vl| with a summary and some speculations 
on the lattice problem. To alleviate cross-referencing, some frequently used parameters and 
symbols are gathered in Table |. 

II. REDUCTION TO ONE DIMENSION AND BOSONIZATION 

In this section, we shall start from the most general two-impurity Kondo model and 
perform various exact transformations to reduce it to a simplified form, fl36|) + fl37|) + fl39|) , 
suitable for identifying the critical point and uncovering the underlying physics. 

The general Hamiltonian for the anisotropic two-impurity Kondo model has the following 
form, 

H = J d*k erffa + J £ £ [e^-^y^ + e-^'y%a^ %l Si\ 

\ / \=x,y,z 

+ K X S*S*, (1) 

\=x,y,z 

where m. = (V'L, wt, ), °" A with A = x,y, z are the Pauli matrices, S\ and S2 are the two im- 
purity spin 1/2 operators located at ±R/2. It has been shown that the Hamiltonian (|TJ) can 



be reduced to an equivalent one-dimensional (ID) problem |[lT| , |18| . Introducing ID fermionic 
operators, 
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with the notations d 3 k = k 2 dk d 2 Q and 
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sin(kR) 



kR 



(3) 



we can completely rewrite Kondo interactions in terms of these new operators 
v F 



H 



Kondo 



2 , ^ _ Jo 2tt 



°° {4(k, k') [4(k)a x Mk') + 4(k)a x Mk') 



Si 



E 

\=x,y,x 

+Ji(k,k') \4(k)a x Mk') ~ 4(k)a x ^ 2 (k')} S x + J x (k,k') \4(k)a x ij 2 (k') 



+4(k)a x Mk')} S x - J x (k, k') [4(k)a x ip 2 (k') - 4{k)a x ^{k')] S x ) . (4) 



In the last expression, we have introduced the short-hand notation, 

S X =S X ±S X , \ = x,y,z. 

The momentum dependent coupling constants are, for X — x,y, z, 

J x kk! 



(5) 
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J x {k, k') 

J m (k, k ) 
J x r (k, k ) 



J x kk' 
4ixvp 
J x kk' 
4i\vp 
J x kk' 
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[N e (k)N e (k') + N (k)N (k')} 
[N e (k)N e (k') - N (k)N (k')} 
[N e (k)N (k') + N (k)N e (k')] 
[N e (k)N (k') - N (k)N e {k')) 



(6) 
(7) 
(8) 
(9) 



From the commutation relation , ipp a ,} = 5 3 (k — k r )5o-y, we can verify 

{4 a (k),^(k')} = 2n5{k - k')5ij>5 a y. 



(10) 



The free conduction electron Hamiltonian can also be written in terms of these ID fermions, 
plus completely decoupled extra degrees of freedom, 

r°° dk 



d 3 k e k ipy 



2vr 



^l(k)Mk) + 4(k)Mk) 



+ 



(11) 



Thus, only ID fermions defined in (0) are relevant to the two-impurity Kondo problem. 

So far, the reduction has been exact. In the next step, we linearize the dispersion at 
k = kp-i €k — vp(k — kp), and expand the momentum dependent coupling constants around 



k = kF- We only need to retain the leading terms of the expansion since other terms contain 
some power of k — kp, which have high scaling dimension and are irrelevant at low energy. 
From J x (kp,kF) = 0, we see that Jj r interaction only contains irrelevant terms. Denoting 
k — kF by k again and with implicit understanding of an ultraviolet cutoff, we can introduce 
Fourier transformations 



30 rjh 

-e^ t (k), z = 1,2. 

oo Z7T 



(12) 



The fermion operators satisfy the standard commutation relation, 



(13) 



After linearization, the full two-impurity Kondo Hamiltonian can be cast in the following 
form, 



H = H + Hi , 

/oo 
dxi/jj(x)d x ipi(x) + K\SiS; 
. , "OO , . 



(14) 



A cA 



i=l,2 ' 



+ K 



1 r°° x 
S z + + 7} E / dx4(x)a*Mx) 

£ : 1 n J —OO 



\=x,y,z 

+ h s S z _, 



i=l,2 



(15) 



^i = y E { J + [^(oyVi(o) +4(0)^2(0) 



X=x,y,z 



+ Ji [^(O)trVi(O) - ^(0)aV 2 (0)] S x _ + J x _ [^(0)^(0) + 4(0)aVi(0)] S\) . (16) 



The coupling constants are, for A = x,y,z, 



J + = IT J p F , J_ = 7T J p F : • 

k F R 



J m = nJ p F * 



( sin k F FC 
\T V kpR , 



(17) 



with p F = k F /(27i 2 vp), denoting the conduction electron density of states per spin at the 
Fermi energy. Noting that we have included both uniform and staggered magnetic fields h u 
and h s in fllBD , with Bohr magneton and gyromagnetic ratio set equal to one. 

At this stage, the 3D two-impurity Kondo model has been successfully reduced to an 
equivalent ID problem, up to some terms irrelevant at low energy. However, we must remove 
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two accidental features of fll4]) resulting from linearization. They are the particle-hole sym- 
metry and a special relation between J±(kp, hp) and J^(kp, kp): that the RKKY interaction 
generated from them is always ferromagnetic ||11| . These accidental features will be spoiled 
by the generated corrections from irrelevant terms neglected during the linearization. That 
an irrelevant interaction can renormalize the coupling constant of a relevant interaction is 
a well known fact p2 |. An example can also be found in section [IV], expression (p2"|) , where 



the dimension 3/2 leading irrelevant operator in the effective Hamiltonian (fyj) induces a 
correction to the dimension 1/2 relevant operator. Usually, the accidental features at the 
lowest order will not survive if there is no hidden symmetry ensuring them. The accidental 
relation between J±(k F ,kp) and J^(k F , k F ), with A = x,y,z, is removed by treating these 
coupling constants as independent parameters. This is also physically meaningful since these 
interactions are completely independent and presumably play different roles at low energy. 
For a general conduction band, particle-hole symmetry breaking, although weak, is always 
present. The general particle-hole symmetry breaking term that can be added to the ID 
Hamiltonian (]TJj) has the following form JTT].[T^ , 



H, = V 



^1(0)^(0) +Vl(0)Vi(0) , (18) 



where V is the energy scale characterizing the strength of particle-hole symmetry break- 



ing. Adding the marginal operator ( lq) to the ID Hamiltonian after dropping irrelevant 
interactions in the linearization may seem unusual, actually it is the natural thing to do. 
The reason is again the generation of (]H|) from irrelevant interactions in the absence of 
the particle-hole symmetry. Usually, all possible operators allowed by the symmetry will be 
generated by irrelevant interactions, and we only need to include relatively more relevant 
operators. In this case, the only non-irrelevant operator breaking particle-hole symmetry 
is ([[8]). In section [TTT] and appendix [FJ, we shall see a similar example where the marginal 
operator (0) generates a relevant operator ( ^7|) when projecting to a subspace relevant for 
the solution region of Figure [TJ 

In the rest of this paper, we shall retain the rotational symmetry around z-axis by setting 

9 



K x = K y = K± and Jf = Jf = for i = m,±. Apart from the continuous U(l) rotational 
symmetry, the Hamiltonian fll^j ) possesses several discrete symmetries which will be very 
useful for our analysis. The transformation rules are, omitting unaffected operators, 



Parity : tpi <-> ip 2 , «-> S$ for A = x, y, z, 
Particle — hole : ip^ — > ipt, ip^ — > —ipjp 
7T rotation around x— axis : ip^ t/}^, Sf — > —Sf, S" — > —S*. 



(19) 



The particle-hole symmetry exists when V = 0. 

The next step is to reduce the Hamiltonian (|14|) + (|18|) to a simple form suitable for 
identifying the critical point. The reduction involves bosonizing the Hamiltonian which only 
contains ID left-moving fermions |23 > p4 , 17 . There are four species of fermions, so we need 
to introduce four bose fields, 



Pr 



where a is the lattice spacing and 



* = 1,2, cr=ti, 



(20) 



^ia(x) = Vtt (f) irT (x) - ( dx' Ii ia (x') 

l J — oo 



(21) 



The bose fields satisfy the standard commutation relation, 



ja (x), U f(T/ (x')} = i 8 jj f5 aa/ 5(x - x). 



(22) 



The phase factors P^ are introduced to take care of the anticommutation relations between 
different species of fermions. Our choices are, 



P 1T = P n = < 
P21 = P21 = ' 



i7T I-00 dx ^lt( x )^lt( x ) 



(23) 
(24) 



By substituting (p0|) into ( |I~4"D and using the relation ipl a (x)ipia(%) = d x § i(T (x)/(2iT), the 
two-impurity Kondo Hamiltonian is expressed in terms of four bose fields <fii a (x). Then we 
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make linear transformations to four new bose fields corresponding to charge, spin, flavor and 
spin-flavor degrees of freedom, 



sf 



hi + <Pn + <M + 02i)/2, 

'IT - 0U + <h\ - 02|)/2, 
'IT + <Pu ~ 02T - 02i)/2, 
'IT - <f>U ~ <M + 02t)/2. 



(25) 



The Hamiltonian now acquires the following form 



#1 



+ K 
v F \ J. 



E r dx {Ul(x) + \dAxix)?} + E K ^i S : 

- - r - r J — OO 1 J , 



2 i ./ — ~\ 

X=c,s,f,sf 



X=x,y,z 



/OO rl-Tf 
-oo Z7T 



2 



(26) 



7* 

v ,, , 



.2 J! 



2Ji 



} ^$ s (0)5; + -J«^$ s/ (0)5! - z— e^cosS^O) sin $/(0)^ + 

2 7T 7T 7T« 7T« 



x cos$ s/ (0) [cos$ s (0)^-sin$ s (0)^ 
2J: 



2JJ 



7TQ; 



sin $ s/ (0) sin $ S (0)S! + cos $ s (0)£^ 



7T« 

.2V 



e l7re sin $/(0) sin $ S (0)S* + cos $ s (0)^ 



# 2 = -i — e my sin$ s/ (0) cos$/(0), 
7TQ; 



(27) 
(28) 



where the phase factor is 
H i .,..r.,.t 



da#L (x)^ u (x) + ^ T (x)^2T(a;)] = 77- / dx[d x $ e (x) - d x $ af (x)]. (29) 

Z7T J-oo 



The charge bose field (j) c {x) is decoupled from the interaction (|27j). It will be omitted from 
now on, and so will be the d x <$> c term inside the integral of the phase 9. The cos$ s (0) and 
sin$ s (0) factors in ( p7j ) can be eliminated by rotating the impurity spins around z-axis by 
an angle $ s (0), 



H -> UHU~\ with ?7 = e - iS +^ {0 \ 

UH,U- l = V -^ J2 r dx{ul(x) + [d x M^)} 2 }+ E K x S^ + h s Sl 



^ C cz ^2 



(30) 



(31) 
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? , ( jz jz 9 J z 

UH.U- 1 = -^$,(0)5; + -^$ s/ (0)S* -i— e*^cos$ a/ (0) sin (0)5; 

2 [ 7T 7T 7TO; 

- -^(S+) 2 + — [J^cos$ s/ (0)5; - J^sm$ sf (0)S y _ - J±e in9 sm$ f (0)Sl] . (32) 

The particle-hole symmetry breaking term H 2 is not affected by the above rotation. We 
note that both J? and K z acquired corrections under the rotation, and h u Sl is canceled out 
in ([H]). The interactions in fl3"2] ) only contain d x <& s (0). Therefore the bose field (f> s (x) can 
be integrated out analytically upon our wish. 

The bosonized Hamiltonian (|31])+ ([32l ) can be re-fermionized by introducing three species 
of fermions, 

iM*) = -l=e^'\ (33) 
V 2ira 

^ f ( x )e i7r f-o°^>^ {x) = — =e^ (a;) , (34) 

V27ra 

4(x)^ s (x) = ±-d x $ s (x). (35) 

Again, a phase factor is included in the definition of the fermion operator ipf{ x ) to take care 
of the anticommutation relations between three different species of fermions. Because the 
interactions in (132]) contain only d x <& s (0)(= 27n/\[(0)?/> s (0)), so we do not need to specify the 



phase for the fermion operator ^ s (0). The complete Hamiltonian can be rewritten as 

/DO . 
dx4(x)d x ip i (x)+K z S z 1 S% + K ± J2 S$S* 

/oo 
dx^l(x)^ s (x) + h s S z _, (36) 
-oo 



Hi = -f { j;[^(o)^(o) - ^(o)^t(o)]s; + ^[^(0)^/(0) - ^(o)^(o)]si 
+j*[tM°) + vI/(o)][^/(o) - mo)]s z + \ + UtlMo) + vI/(o)]5; 

+iJi[^/(0) - V>I/(0)]# + iJ±#/(0) - ^ f (0)]Sl} , (37) 



where 



j; = j; - 2vr, K Z = K Z - — (j; - tt) . (38) 

7TQ? V 7 



The particle-hole symmetry breaking term becomes 
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Ho = V 



V> s/ (o) - V>l/(o)] k/(o)+^(o) 



(39) 



How do these new fermion operators transform under the discrete symmetries of (|T9|)? 
We can keep track of the transformation rules during the bosonization and subsequent 
fermionization to derive, omitting unaffected operators, 



Parity: ip sf ipl f , $ f 
Particle — hole : ipj *-> —ipt, 



S$ <-> for X = x,y, 



(40) 



7r rotation, x— axis : ip s f ip\ f , ip s <-> Vl> ^ ~~ — > S? — > — S£ 



Alternatively, one can directly verify them from (|36|) , (|37|) and (|39"|). 



III. EFFECTIVE HAMILTON! AN NEAR THE CRITICAL POINT 

In this section, we shall identify the critical point from (p6|) + (p7|) an d derive an effective 
Hamiltonian for the finite region of the parameter space surrounding the critical point. Our 
identification of the critical point will make clear its physical origin. 

To search for the critical point, we need only consider the particle-hole symmetric case, 
V = 0. At first glance the Hamiltonian (|36| ) + (|37|) still looks too complicated to provide any 
intuition. On the other hand, from the conformal field theory results it is known that the 
critical point exists in a restricted Hamiltonian with Jl — J_ — 0, and around the critical 
point J*, J_ interactions are irrelevant. Thus, our task is greatly reduced by searching the 
critical point in this restricted Hamiltonian. An important step is to verify the irrelevance 
of the Jf, J^z interactions after we find the critical point. 

We have noted before hat only the product ipl(0)ip s (0) appears in (|37|). This is because 
only d x <& s (Q) appears in ([32]) and there isn't anything containing cos $ s (0) or sin $ s (0). Thus, 
the Jf term containing the bose field <f> a (x) can be integrated out analytically. Although 
d x <& s (0) couples to an operator SI., the integration can be done formally in the path integral 
formalism, yielding the following two terms to the action, 
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- V -4^r- f dr[Sl{r)] 2 + i^±?J2Wn\S z + (-v n )S* + (u n ), with u n = 2mr/f3. (41) 

The first term is a correction to the RKKY interaction, so it is absorbed into K z . The second 
term has higher dimension and is expected to be irrelevant. The point we want to make 
here is that the J\ interaction does not affect the critical point and can be ignored during 
the search for the critical point. Thus, we see that when J z _ = J^; = the remaining Kondo 
interactions only involve three local spin operators, S+, SZ, St, which only act on three (out 
of four) local impurity spin states. The impurity spin state (| |T> — I ll>)/v / 2 decouples 
from the Kondo interactions when J z _ = J_ = 0. Together with the RKKY interactions, we 
derive an energy level scheme for the impurity spin states in Figure The critical point 
corresponds to the special case when the two lowest levels become degenerate. Specifically, 
the —K z (S : L) 2 /2 term in (|36|) raises the energy of the states | f|> and | |f> by an amount 
—K z /2 (assuming — K z > 0) with respect to the other two states | ft> ±| ll>- These two 
states, | tl> an d | IT>> are further split symmetrically by the transverse part of the RKKY 
interaction, K L (StS 2 + S{S£)/2. When -K z = K ±} the two levels, (| ||> - | IT>)/v / 2 
and (| ft> + I ll>)/v / 2, become degenerate, forming a doublet. Because there is almost 
no Kondo interaction in the state (| |T> _ I l|>)/v / 2 when J z _ = J 1 : = 0, the superficial 
degeneracy between this state and the doublet is lifted by the Kondo interactions in the 
doublet which lower the energy of the doublet by a finite amount T K , equal to the ground 
state energy gain at the critical point. Turning on J z _ and J_ will not change the energy 
level scheme as long as vpJ^, vpJ^; < T K . 

Now, we have identified the critical point. To describe the low energy physics, it is 
sufficient to project the full interacting Hamiltonian (|36|) + (|37|) + (|39|) onto the lowest energy 



doublet. For J^_, J^, J^, J* , J^, V all but much smaller than one, the projection can be 
done accurately. Let Q be the projection operator onto the doublet, the projected effective 
Hamiltonian to the second order is 

g„, = QHQ + QH (1 - Q) Eo . RQ 1 ^ q) (1 - Q)HQ, (42, 
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where E$ is ground state energy. The doublet can be described by local fermion operators 
d and d) such that (| ft> + I ||>)/v / 2 and (| ||> - | IT>)/v / 2 correspond to d ] d = and 
d)d = 1 states respectively. From Figure |3|, it is not difficult to verify QS^LQ = i(d — eft), 
QS\Q = QS+Q = QS±Q = 0. These relations are used for the first order projection. In 
the second order projection, nonvanishing terms may contain Q(Si) 2 Q = d)d, Q(S rJ l) 2 Q = 
Q(S^) 2 Q = dd\ QS+S Z _Q = d and QS^S^Q = d^. Since the extended fermions commute 
with the impurity spin operators, we need to to install anticommutation relations between 
the local fermion operators d, S and the extended fermion operators ip\(x) with A = s, f, sf. 
This is accomplished by a simple transformation, 



ip\{x) = ip\{x) e 1 



TdU 



A = s, f, sf. 



(43) 



The commutation relations between d and ip\{x) are converted to anticommutation relations 
between d and ip\(x). The effective Hamiltonian will be represented in terms of d and 
ijj\(x). But we shall omit the tilde signs on ip\{x) in the following. With the help of the 
above mentioned results, it is straightforward to evaluate (0). The details are captured in 
appendix O. The results are 



Heff 



H, 



-iv F 

i=f,sf J - co 

+a u h u [Vv(0) + <(0)] ty/(0) - 4(0) 
K ± + K : 



pert 



Hfp + Hp er t + H p hb, 

coo r -, 

dxipl(x)d x tf; i (x) + v F g ip sf (0) - 4/(0) J (d + d ] ) 

+ a s h s [^ s /(0) + 4/(0)] (d-d^), 
K c ) d ] d - ivpgx (d - d ] )d x ip s f(x) - i>l f (x) 



x=0 



^ = ^[^/(0)+4(0)](rf-rf t ), 



(44) 

(45) 
(46) 
(47) 



where K c is the critical value of [K z + K±)/2. We have separated the Hamiltonian into the 
fixed point part, perturbation part, and particle- hole symmetry breaking part. We note that 
(pE5j) in the absence of the external magnetic fields h u and h s has the same form as that of 
the two-channel Kondo model f24|. To the second order, the coefficients in H e ff are given 
by 
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Ji 

go = 



V2 



ira 



in 



4Tra(K ± +T K )J^ 



(48) 



„2 7_L jz 

0l = V F-J+J m (4Q) 

yi 8(ir ± + T^) 2 (27ra) 3 / 2 ' 1 ' 

V = V F J + J m V /roN 

2(K ± + T^) 2 (2tt«) 5 / 2 ' 1 J 

°« = "g^' (51) 

= W + - (52) 

(K ± +T K ) v / 2^ 

The projection induces corrections to the RKKY interactions so the critical value K c is 
not exactly vf(J+ — 7r)/(7ra), as determined by the condition — K z = K±. The energy 
Tk in the above expressions can be identified as the ground state energy gain in (|45|) from 
the hybridization term Vpg$. From the study of the two-channel Kondo problem, we know 
T K ~ v F gl. 

The spin degrees of freedom are completely decoupled and their Hamiltonian is, parallel 
to O, 

/oo poo Vt?( J z \^h 

dx^\{x)d x ijj s {x) + h u / dxiP\{x)^ s {x) + F \ % u 4 (0)^(0). (53) 
-oo J-oo 47Ti^ 

Since this piece of Hamiltonian does not contain any interesting physics, we shall not discuss 
it hereafter. For the particle-hole symmetry breaking and staggered field coupling terms, 



apart from the relevant operators we have included in (f45|) and (f|7|), there are also marginal 
operators such as d39|). They will be discussed below and in appendix [TJ. 

The effective Hamiltonian ( f44|) is the central result of this paper. The rest of this section is 
to prove the completeness of the (|44] ) from general symmetry considerations. For a certain 
region of the parameter space, the projection is controllable in a sense that high order 
corrections to the coefficients of H e ff are too small to alter its critical behavior. However, 
the projection is not done exactly. One may ask how do we know that there are no other 
operators which could arise from high orders of the projection and spoil the critical behavior? 
Fortunately, it turns out that all other operators up to dimension 3/2 inclusive can be 



eliminated by the three discrete symmetries of (^0|) . To show this, we first determine how 
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the operators d and d) transform under parity and the rotation of an angle 7r around x-axis. 
The particle-hole transformation does not involve the impurity spins, so will not affect d 
and d\ Since the two states of the doublet have different parity and d\ d connect them, we 
conclude that under parity: d) — > —d\ d — ► —d. This could also be seen from QStS+Q = d) 
and QSIS^Q = d. From them, we can also see that under the ir rotation around x-axis, 
d) — > —eft, d — > —d. Combining these results with (f5]), we derive all "elementary operators" 
in the projected Hilbert space comprised of the local doublet and the extended fermions 
ip\(x) with A = s, f, sf. They are listed in Table |T|. 

Some explanations are necessary at this point. First, it is well known that ID extended 
fermion operators have scaling dimension 1/2. This can be easily seen from the free fermion 
action S(ip, ift) = Jq dr f^° 00 dx^(x,r)(d T — ivjrd x )ip(x,r). Secondly, we have noted before 
that ip s and ip\ can only appear in the product ip\ ip s . Therefore, the spin degrees of freedom 
do not bring in the dimension 1/2 operators ip s and ip\ as additional building blocks in 
Table 0. Thirdly, usual local fermion operators have scaling dimension zero. As can be seen 
from the free fermion action S{d, (ft) = fft drd) (r)d T d(r), we need not change d and d) under 
a rescaling of the imaginary time r. This would imply that both combinations d ± dJ have 
scaling dimension zero. However, the dimension of the operator combination d + d' is raised 
to 1/2 by the hybridization term in fl45|) with coefficient vf9o- This follows immediately from 
the requirement of preserving scale invariance of the hybridization term under a rescaling 
of x and r. Lastly, because both d and d) are odd under parity and 7r rotation around 
x-axis, only ipsf{0) — ^1/(0) could hybridize with them to give rise to a term in the effective 
Hamiltonian even under all discrete symmetries, as can be seen from Table [n[ Out of two 
linear independent combinations from d and d) , i/j s f(0) — ^1/(0) could only hybridize with 
one. Requiring it to be hermitian, the hybridizing combination could be either d + gF or 
d - d). Thus, once [ip sf {0) - ^ sf (0)](d + d^) is generated in (g5j), (d - cft)[ip af (0) - ^(0)] is 
forbidden. This guarantees that the dimension of d — (ft will remain to be zero, like usual 
local operators. 

To construct operators in the projected Hilbert space, we only need to multiply together 
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the building blocks of Table and keep products of even number of fermionic operators. 
We list all dimension 1/2 operators in Table [HT] , all dimension 1 operators in Table [TV], and 
all dimension 3/2 operators in Table |V|. 

Let us first consider the particle-hole symmetric case. Any operator that could appear in 
the effective Hamiltonian must be even under all three discrete symmetry operations (flOp. 
We can explicitly verify that all allowed operators up to dimension 3/2 that could appear 
in the effective Hamiltonian have been included in (flq ) and (|46|). In order to couple to the 
uniform magnetic field, an operator has to be even under parity and particle-hole transfor- 
mations but odd under 7r rotation around x-axis. From Tables |T| and [TV|, we also verify 
that the only allowed operator up to dimension 1 inclusive is the one appearing in (|45|). As 
to the operators that could couple to the staggered magnetic field, they must be even under 
particle-hole transformation but odd under parity transformation and 7r rotation around x- 
axis. Apart from the dimension 1/2 operator that couples to the staggered field in ([45]), two 
more dimension 1 operators are allowed by the symmetries. They are the third and fourth 
operators in Table [TV]. Thus, including marginal operators we could have the following 



additional staggered field coupling terms added to (01) 



Hi 



stag 



a% k/(0) + 4 f (0)] k/(0) - 4 f (0)] + ia" 8 h s k/(0) + ^(0)1 (d + tf), (54) 



where a' s and a" are two dimensionless parameters depending on the original coupling con- 
stants of QT4]). Nevertheless, the contributions to the staggered susceptibility from (|54]) are 
negligible as we shall see in appendix |D|. 

A subtle point arises here. We have used the argument, that Vv(0) — ^1/(0) could 
hybridize with either one of d ± S but not both, to rule out possible hybridization between 
■0s/(O) — ^1/(0) and d — d^ in (^). Why couldn't we make the same argument to eliminate 
[ip sf (0) + ipl f (0)](d + d^) in (H, since we already have [V>./(0) +V>J / (0)](d-d t ) in (||). The 
reason is the following. For an arbitrary hybridization, we can rewrite it in the following 
way 



^ s /(0) - -01/(0)] [a(d - + i/3(d + S)] =^/ a 2 + f32 [Vv(0)-V>!/(0)] (de^-Se 



-tip 



where a, j3 are two arbitrary real constants, ip = tan _1 (/3/a), and the % in front of j3 in the 
left side of the last formula is needed to make that term hermitian. Redefining the operators 
d and gF to absorb the phase (p, we reduce the hybridizing combination to either d — eft or 
d+d>. In other words, we can always choose a proper definition for d and d' so that only one 
of d±cF hybridizes with ^ s /(0) — 4/(0)- But we can only perform phase absorption once. A 
redefinition of d and d) to absorb a second phase to eliminate [ip s f(0) + 4/(0)](e£ + o^) from 
the stageered field coupling terms in (^5|) + (|54|) is not possible without spoiling the simple 
hybridization form of the fixed point Hamiltonian (|45f) . 

When the particle-hole symmetry is broken, another relevant operator becomes allowed, 
as can be seen from Table [II. This is the dimension 1/2 operator in (|47|). There are also 



two dimension 1 operators breaking only particle-hole symmetry. From Table |V|, they are 

Km = v K/(o) - 4/(o)] [v>/(o) + 4(o)] + z^v [V/(o) + 4(o)] (d + t/t), (55) 

where a v is a dimensionless coefficient depending on the original coupling constants of ( |14| ) . 
The first term in fl55|) is the original particle-hole symmetry breaking term ([39|), surviving 
the first order projection. The second term is a generated one from high orders and cannot 
be eliminated by a simple phase absorption in d and d^ for the same reason of the last 
paragraph. 

Summarizing this section, (f4l|) + fl54|) + (|55D constitutes the most general effective Hamil- 
tonian for the solution region of Figure |I|, even allowing particle-hole symmetry breaking. 
What are omitted up to dimension 3/2 only include: 



A dimension 1 operator Vv(0) +4/(0) ^/(O) +4(0) > which could couple to the 
staggered field h s but breaks the particle-hole symmetry. Thus, the coefficient of this 
operator must be proportional to the particle-hole symmetry breaking potential V. 
Close to the critical point, we expect this coefficient to be significantly suppressed. 
This term should be even less important than those in (0). 

Two dimension 3/2 operators breaking only particle- hole symmetry, as can be seen 
from Table 0. They could appear as additional irrelevant operators in the effective 
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Hamiltonian. Again, we expect they are significantly suppressed close to the critical 
point. 

• Several dimension 3/2 operators which could couple to the uniform or staggered mag- 
netic fields. Their contributions to the susceptibilities vanish according to high powers 
of temperature as T — > 0. 

It is worth pointing out that up to dimension 3/2 the number of allowed operators around 
the critical point and their dimensions are in complete agreement with the conformal theory 
results |i~8|| . 



IV. LOW ENERGY THERMODYNAMICS 

In this section, we shall calculate low energy thermodynamic properties of the effective 
Hamiltonian for the solution region of Figure [I]. The marginal operators ( pi] ) will be 
considered in appendix 0, where we shall show that their effect is to slightly renormalize the 
Kondo and crossover temperatures. The contribution to the staggered susceptibility from 
the marginal operators ( 54j) will be considered in appendix [D] and shown to be negligible. 



The way we shall adopt to carry out calculations is to represent the partition function as 
a path integral in which every fermion operator becomes a Grassmann variable. Then we 
perform linear transformations on the Grassmann variables to bring the action to a diagonal 
form. 

The partition function in the path integral formalism can be written as 

Z = J V[i; sf ^ sf ,^ f ,^ f ,d,d]e~Io d ^ +H ^\ (56) 

/■oo r _ _ - _ 

C = dx ip s f(x)d T tp s f(x) + ipf(x)d T tfjf(x) + dd T d, (57) 



where H e ff is the effective Hamiltonian given by (0). By making linear transformations to 
new Grassmann variables, 



1 % 

-j=(i) sf + i) sf ), & s/ = -— v 



hf = —fi=(lf>af + ll>3f), b af = 7=008/ - VVf)) ( 58 ) 
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a f = + 



V2 



--^{ipf-tpf), 

(d — d), 



we write the total Lagrangian in fl56|) as 



(59) 
(60) 



C = C + H eff = £i(a s/ ) + C 2 (b sf ) + C 3 (a f ) + C A (b f ) + Ci oc {a, b), 

] roo 



£i( a sf) = o / dx a sf (r,x)(d T - iv F d x )a sf (r, x) + 2ia s h s a s f(r, 0)fe(r) 

2 J—oo 



C 2 (b sf ) = - dx b sf (r,x)(d T - iv F d x )b sf (r,x) 

2 J—oo 

+2iv F [g b sf (r, 0)a(r) + gib(r)d x b sf (T, 0)] , 



jC 3 (aj) = — / dx a/(r, x){d r — iv F d x )af(T, x) + 2iVaf(r, 0)6(r), 

2 J—oo 

1 r 00 

= 77 / dxbf(T,x)(d T -iv F d x )bf(T,x) + 2ia u h u a s f(T,0)bf(T,0), 

2 J—oo 

C loc (a, b) = i [a(<r)cU(-r) + fc^)^)] + * ^ a(r)6(r), 



(61) 
(62) 

(63) 
(64) 

(65) 

(66) 



where 



air = # c + id). 



Introducing the Fourier transformation, 



./( 



a,HT,x 



-iu} n r+ikx 



(67) 



(68) 



and similar transformations for the other Grassmann variables, we can write the actions 
corresponding to the Lagrangians fl62|) to flBTf ) as 

1 



Si{a sf ) = E/ 



00 dk 
oo 2tt 

+2ia s h s a sf (uj n , k)b(-u n )] 
00 dk 



iu n - v F k)a sf (-uj ni —k)a s f(uj n , k) 



/OO (Jfe f I 
— <^ --(iu n - v F k)b sf (-u n , -k)b sf (u n , k) 
-co 2lX I 2 

-2iv F [g Q a(-uj n ) - igt k b(-w n )] b sf (uj n , k)} , 



dk 



oo 

-oo 2n 

°° dk f 1 



--(iuj n - v F k)a f (-ui n , -k)a f {u n , k) + 2iVa f (-u ni -k)b{u n ) 



n 

/oo r | 

— -~(iu) n -v F k)b f (-uj n ,-k)bf(uj n ,k) 
-oo 2ix L I 



(69) 

(70) 
(71) 
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+2ia u h u a sf (-u n , -k)b f (uj n , k)} , (72) 
Si oc (a,b) = [a(-uj n )a(uj n ) + b(-u n )b(u n )] +i5K a(-u n )b(u} n )j . (73) 

The uniform magnetic field term in (|72|) is an exactly marginal operator. This is most easily 
seen by setting h s = and combining ( |6~9|) with (|72"|) . From the real Grassmann variables 
a s f and bf, we can make a linear transformation to 



^ = 71 + ^ = 71 ~ 



The Grassmann variables ?/> and ^ correspond to the usual fermion annihilation and creation 
operators. When h s = 0, (|69|)+([7^) is completely decoupled from the rest of the Hamiltonian 
responsible for the critical behavior. Moreover, the h u term of ( |72|) is simply a potential 
scattering term in terms of the fermions corresponding to ip and 



/■p ("OO _ _ 

S(a sf ) + S(b f ) = dr dxip(x)(d T - iv F d x )ip(x) + 2a u h u ^ (0)^(0) 

JO U-oo 



(74) 



Thus, not only the uniform susceptibility is well behaved but also applying a uniform external 
magnetic field has negligible effect on the physical behavior of the system inside the solution 
region of Figure [l|! From now on, we shall set h u = and drop ([72]) from further discussion. 

We can diagonalize ( |69| ) to (|7T|) simply by shifting the Grassmann variables corresponding 
to the extended degrees of freedom, 

a sf (uj n , k) = a sf (uj n , k) - 2 ^A ^ ) (75) 

6 s /(o; n , k) = b sf (u n , k) - . 2Wp - [g a(u n ) + igxkb^n)] , (76) 

iuj n — v F k 

2iV 

a f {u n ,k) = a f (uj n ,k) - - rb(u n ). (77) 

iuj n — v F k 



Upon inserting the results for the following integrals, 

dk 1 i sgna; n 

-oo 2tt iuj n — v F k 2vf ' 

A dk k 1 / vpA lav, 



a 27T iui n — Vpk V F \7T 2 



J F 

„2 



A dk k %uj n I VpA \uj rt 

—A 2tt iuj n — Vpk vp \ it 2 
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(78) 
(79) 
(80) 



where A is the ultraviolet cutoff, the actions (|69D to (1711) become 



Si (a. 



sf) 



E 



00 dA; 

-00 47T 

i{<y s h s ) 2 

v F 
r °° dk 



[iu n - v F k)a sf (-uj n , -k)a s f(u n} k) 
^2sgnu n b{-u) n )b{u n ), 



(81) 



— (ioj n - v F k)b sf (-uj n , -k)b sf (uj ri , k) - iv F gl Y sgnw n a(-a; n )a(a> n ) 

-00 47T „ 



+2ig gi Y] \u n \a(-uj n )b(uj n ) + -Vw„ (|a> n | - ) b{-u n )b{u n ) 
4iv F g giA 



7T 



J^a(-u; n )6(u; n ). 



(82) 



5 3 (a/) 



^ f°° dk. 



iV 2 



v F k)a f (-uj n , -k)af(u n , k) ^ sgnu;„ 6(-a; n )&(a> n ). (83) 



The last term in (|82 ) is a correction to the relevant operator term of (|73|) and can be 
absorbed into the critical value of the RKKY interaction K c . Collecting the local terms 



containing Grassmann variables a, b from (|5ID to (|83| ) and combining them with (fT3|), we 
obtain the effective local action, 



s io f I = - l E {a(-uj n ),b(-uj n )) 



n>0 



X 



\u n \ + 2v F gl 
5K + 2g g 1 \u n \ 
The factor is defined by 



V 



\0J, 



-(8K + 2g gi\uj n \) 
„\/Z b + 2(a 2 s h* + V2-g!u;l) /v F ) 



b(u r , 



(84) 



(85) 



1 + 4^ 2 A/7t' 

and can be interpreted as the wave function renormalization factor for the Grassmann vari- 
able 6(or Majorana fermion). The a; 2 term in the matrix element of (|S4|) can be safely 
neglected since it is highly irrelevant and satisfies giuj^/vp <C \oJ n \ for the whole energy 
range of practical interest. All interesting thermodynamics is contained in ([S4]). From ( |84|) 
we obtain the free energy shift due to the local interactions, 



F(T, h s ) = -Tln2-Tj2 hi{{6K + 2g ogi \u; n \f 



n>0 



+ (K|+2^ 2 ) ^ + 1(^ + ^)1} 



(86) 
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where — T In 2 is the entropy of two degenerate impurity spin states. Defining several con- 
venience notations, 

T K = 2v F gl + 2Z b + 2g o9l 6K^j , (87) 

a s = 2g a s JZ b , (89) 



we can recast the free energy in a very simple form, 

V ' ; 7-do 2tt e#" + 1 \ cj 2 - T C T K - a 2 s h 2 s J v ; 



The roles of the parameters in (9~0) can be read off. Tk is the fundamental energy scale of 
the problem and should be identified as the Kondo temperature. We note Tk — 2vf9q, as 
can be seen from ( j37|) . T c vanishes approaching the critical point, and satisfies T c Tk 
inside the solution region of Figure |I|. The same T c defines the crossover energy scale above 
which the behavior of the system is controlled by the critical point. Below T c , it is controlled 
by the Fermi-liquid fixed point. Accompanying the staggered magnetic field is an involved 
coefficient a s because h s couples to an unconserved operator. Because of this factor, it is 
not possible to define a universal Wilson ratio from the staggered susceptibility. 

For all practical purposes, the h 2 s term in the numerator inside tan -1 in ( p0|) can be 
dropped since it only shifts Tf to T K + a 2 s h 2 s . After some rearrangement, we finally obtain 



This nice looking expression gives us the complete crossover functions for the specific heat 
and staggered susceptibility. 



A. Specific Heat 

Setting h s = in ( |9~l"D and performing some minor manipulation, we obtain 
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F(T) — F(0) = -Tln2 + T / 

J o 



00 dx 1 



tan 



-1 



xpTj 



K 



lo 7T e x + 1 V^ 2 - /5 2 ^ C 7V 

Two limiting behaviors follow immediately. At T c <C T <C T^, 

f°° dx 1 



re 

F(T) — F(0) = -T\n2 + T / 

Jo 



71 X 



T Ti T 

— In 2 . 

2 12 



vr e x + 1 \2 /3T K J 
We see a residual entropy (ln2)/2, reduced from the original In 2. At T <^T C <^ T K , 

r°° dx 1 / 1 \ 7r 1 



F(T) -F(0) = -Tln2 + T / 

Jo 



71 — 



ti e x + 1 V" /3T C ; 12 T ( 
The limiting behaviors of the specific heat are obvious from (|93|) and (0 

C(T,T C ) _ f efc, T»T C , 



-T 



T 



— T ~» T 



The general crossover function for the specific heat is obtained from (|92] 



C(T) ^ 2 n4f°°dx x x 4 T K - T c (x 2 - f3 2 T c T K )(3x 2 + f3 2 T c T K ) 



(92) 



(93) 



(94) 



(95) 



(96) 



T """Jo 7re* + l [(x 2 - (3 2 T C T K ) 2 + (x(3T K ) J 
This crossover function is plotted in Figure [|. We also plot the crossover function for the 
entropy ( |9"2"D in Figure [5], which may be more instructive. 



B. Staggered Susceptibility 



From (|9lD, the staggered susceptibility is given by 



Xs (T,T c ) = - 



d 2 



F(T,h s 



h s =0 



oc 



du 



'/V 



00 <T, 



K 



The limiting behaviors are found to be 



Xs(T,T c 



ttT k 



MnT, T>T C 



ttT k 



t-lnT c , T<T C 



(97) 



(98) 



The crossover behavior for the function (|97|) is plotted in Figure || 
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C. Impurity Spin Correlation (S\ ■ S 2 ) 

Projecting the operator (S+) 2 onto the doublet, we find 

QS 2 + Q = Q(S X + ) 2 Q + Q(S y + ) 2 Q + Q(S Z + ) 2 Q = 2dS . (99) 

Thus, the calculation of the impurity spin correlation is reduced to evaluating (deft). From 
( jS4]) , we have 

(ddt) _ I = i{ba) = TE 6K+2go 9l \u; n \ = 

2 ^ (5K + 2g ogi \uj n \y + (\u n \+T K )(\uj n \/Z b + 2Vyv F ) 

The impurity spin correlation is 

<&•&>= 4 +#s "^''"i • HOD 



At the critical point, 5fC = T c = 0, and as T — > 0, 

(3l • 3 2 ) ^ --r + Z b — - -7 + ln ^- ( 102 ) 

The leading irrelevant operator induces a small non-universal correction to the impurity 
spin correlation. Although the fixed point itself (the ground state) has an extra symmetry 
between the two states of the doublet, d <-> d\ which implies (Si ■ S 2 ) = —1/4, it is broken 
by the leading irrelevant operator, as can be seen from (f46|). 

We can also calculate the slope of the impurity spin correlation with respect to the 



RKKY interaction. From (101), we find 



d Z b ^ 1 



d (5K) z/ Pfu* + \u n \T K + T C T K 



00 



M T tanh {^) {Tc T K -$r + (w (103) 



Comparing the last expression with (|97|), we find 





(Si- S 2 ) = const. x Xs (T,T c ). (104) 



d(5K) 

In particular, they should have the same limiting behaviors as given by ( |98|) . Since 

° 0i ■ S 2 ) ~ f dr(f(S x ■ S 2 )(t) 1 ■ S 2 )(0)), (105) 



d (5K) 







the result (|104|) should not be too surprising. 

26 



V. COMPARISON WITH OTHER RESULTS AND UNIVERSALITY 



There are two kinds of asymptotically exact limiting results with which we can compare 
our solution. These are the conformal field theory results at T > T c and the numerical 
renormalization group results at T = 0. First of all, we would like to emphasize that all 
our results are also asymptotically exact up to some numerical coefficients go, g\ and V in 
([H]). Or eventually, the possible uncertainty boils down to the two basic energy scales T# 
and T c , which we cannot determine exactly in terms of the initial parameters of the original 
Hamiltonian ([Tp. 

In order to make a comparison, we first need to answer the questions that whether or not 
the critical point we have studied is the same one, and whether or not the spin anisotropy 
we have introduced in is irrelevant. The answer to both questions is a convincing yes, 
if not rigorous. A detailed comparison of the finite size spectrum of the critical point 
between the conformal field theory and numerical renormalization group approaches has 
been made Excellent agreement has been found which implies the same critical point 
in those two approaches. Thus, we shall take the agreement between our results and that 
obtained from either one of those two approaches as a positive evidence. The conformal 
field theory tells us that there is only one non-Fermi-liquid fixed point, i.e. conformally 
invariant boundary condition [IBf. This is supported by the failure of finding other critical 
points in our approach by considering other impurity spin states as the lowest degenerate 
levels than the doublet of Figure |2j. The strongest evidence for the universality of the 
critical point is the exactly same operator content around the critical point we find in our 
approach and in the conformal field theory approach. That is we have the same number of 
operators with the same symmetry and same dimension. Specifically, there is one dimension 
1/2 relevant operator and one dimension 3/2 leading irrelevant operator in the presence of 
the particle- hole symmetry, as can be seen from (f|5|) and (|46|). Breaking the particle- hole 
symmetry introduces another dimension 1/2 relevant operator, as can be seen from (|47l). 
Furthermore, the dimension 3/2 leading irrelevant operator in the conformal field theory is 
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a descendent of the relevant operator. In our approach, we consistently find that the leading 
irrelevant operator contains the spatial derivative d x . This is the crucial difference from the 
two-channel Kondo problem, resulting in different low temperature behavior for the specific 



heat (see ( fffi) and reference ||24|| ). While in both cases there is a dimension 3/2 leading 
irrelevant operator, only the one at the critical point of the two-impurity Kondo problem 
contains d x . As to the spin anisotropy, it is found in the conformal field theory approach that 
a small spin anisotropic perturbation around the critical point is irrelevant []18| |. Although 



this does not prove the irrelevance of the spin anisotropy introduced in our approach because 
the introduced anisotropy is not small, it does point to the right direction. It is worthwhile 
recalling that the spin anisotropy is irrelevant for all kinds of one-impurity Kondo problem, 
including the exactly screened |23| and overscreened cases ||25| , p6jj . Leaving aside specifics 



of the two-impurity Kondo model, these early experiences in related Kondo problems give 
us considerable confidence that the critical point of the two-impurity model that has been 
studied from different approaches belongs to the same universality class. 

Since the behavior of the system above T c and the way the system flows to the stable 
Fermi-liquid fixed point below T c are all governed by the critical point, the universality of 
the critical point also implies the universal behavior inside the solution region of Figure [I], 
as well as all crossover functions. In particular, the crossover functions we have derived for 
the specific heat and staggered susceptibility are expected to be universal. For a comparison 
of the results from different approaches, the only freedom left is to match the two basic 
energy scales T K and T c . For the staggered susceptibility, or any other response function of 
a non-conserved operator, there may also be an undetermined overall constant prefactor. 

At the critical point T c = 0, or more generally in the limit T ^> T c , the critical properties 
of all thermodynamic quantities as a function of the temperature that we have calculated in 
the last section completely agree with the conformal field theory results, as expected on the 
grounds of the same operator content. These include: the residual entropy (ln2)/2, linear 
specific heat, InT singularity in the staggered susceptibility, constant uniform susceptibility, 
and InT singularity in the correlation function of the composite operator Si ■ S2, as can be 
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seen from ( |103|) and ( |105|) . The complete agreement of the critical behavior further ensures 



us the universality of the critical point. 

As to compare with the numerical renormalization group results at T = 0, we first note 



that the empirical observation [Tl[] of an additional hidden symmetry between the singlet and 
triplet impurity spin states at the critical point becomes crystal clear after our identification 
of the critical point, as can be seen from Figure ||. So is its consequence about the value of 
the impurity spin correlation at the critical point, (Si ■ S2) = —1/4. However, this hidden 
symmetry is broken by the leading irrelevant operator, as we have noted before. This is 
similar to the well known dynamical symmetry breaking effect. In consistency with the 
numerical renormalization group result, we also find that the linear coefficient of the specific 
heat diverges quadratically in 6K on the particle-hole symmetric axis, as can be seen from 
(p5|). Our result ( |104|) , that the slope of the impurity spin correlation with respect to the 
variation of RKKY interaction is logarithmically divergent, is broadly consistent with the 
numerical renormalization group result which also found it singular. 

The only disagreement with the reported numerical renormalization group results is the 
behavior of the staggered susceptibility at T = 0. While Xs(T = 0,T C ) ~ 1/T C has be 
claimed, we only find Xs{T = 0,T C ) ~ hiT c , as can be seen from (|9lf). Noting that the 
other limiting behavior of (|9~8l) , Xs(T,T c ) ~ InT at T ^> T c , is not disputed. Even if 
one takes a cautious view about the numerically fitted critical exponent 2 for the stag- 
gered susceptibility, i.e. Xs ~ ($K)~ 2 , the original numerical divergence seems to us much 
stronger than a logarithmic singularity. The reason for this discrepancy is unknown at this 
moment. But at least the easy explanation of differently adopted definitions for the stag- 
gered susceptibility is unlikely. In this paper, we only couple the staggered field h s to St 
in (|T5|). One could also couple h s to [ip\(0)a z ipi(0) — ^2(0)o" 2 V'2(0)]/2 in (fTJ), or even to 



/ dx[ipi(x)a z ipi(x) — ■ip2( x )°' z ' l p2(x)]/2. In any case, the contributions to the staggered sus- 
ceptibility after subtracting out the free Fermi sea contribution should only come from the 
local operators which are odd under parity and 7r rotation around x-axis. Since we only 



have one such relevant operator, as can be seen from Table |T|, we do not expect qualitative 
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change of the behavior of the staggered susceptibility as a result of different definitions. 
A careful reexamination in the numerical renormalization group approach should be very 
helpful to clarify this point. 



VI. CONCLUSION 

We have presented an asymptotically exact solution for the two-impurity Kondo model 
for a finite region of the parameter space surrounding the critical point, as shown in Figure |l|. 
We have also derived the analytic crossover functions for the specific heat and staggered 
susceptibility. This solution is made possible by an explicit identification of the critical 



point and its underlying physics. As we have explained in section |TJ, the condition for the 
criticality is the degeneracy between the two lowest impurity spin states, (| tl> - | IT>)/\/2 
and (| ft> + | ||>)/v / 2- By varying RKKY interaction across the critical point, these two 
levels cross each other. The non-Fermi-liquid behavior at the critical point is a consequence of 
the fact that in the presence of the particle-hole symmetry only half of the degrees of freedom 
of the doublet can be compensated by the extended degrees of freedom associated with the 
conduction electrons. Because of the special symmetry of the doublet, i.e. one level is an even 
triplet and the other is an odd singlet, the local degrees of freedom of the doublet can only 
be compensated by the extended spin-flavor degrees of freedom of the conduction electrons. 
Out of four species of spinless fermions or eight species of Majorana fermions associated with 
all conduction electron degrees of freedom, only one species of Majorana fermions, i/j s f—ip s f, 
is allowed by the symmetry to compensate the local degrees of the freedom of the doublet. 
This is the same physics responsible for the non-Fermi-liquid behavior of the two-channel 
one-impurity Kondo model. However, the doublet at the critical point of the two-impurity 
problem has different symmetry from the simple impurity spin up and down states of the two- 
channel problem. Therefore, the operator contents around the fixed points(one unstable, the 
other stable) are different. The nearly complete agreement of our results with those derived 
from the numerical renormalization group or conformal field theory approaches, except one 
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limiting behavior of the staggered susceptibility, convincingly establishes the universality 
of the critical point. Thus, the crossover functions we have derived in section [IV] are also 
expected to be universal. The calculation of dynamical correlation functions such as the 
conduction electron Green's function is currently under way. 

What have we learned about the lattice problem from the study of the two-impurity 
Kondo model? An obvious lesson is learned from the striking difference between the uniform 
and staggered susceptibilities. This difference is solely due to the competition between 
RKKY interaction and the Kondo effect. A direct and primitive translation to the Kondo 
lattice problem would be the strong momentum q dependence of the spin susceptibility 
x"(u;, q). As a result of the competition, we should expect drastically different enhancement 
at different momentum transfer q. From this perspective, the picture of a periodic array of 
coherent Kondo scattering centers for the heavy fermion compounds is surely oversimplified. 
Non-perturbatively incorporating RKKY interaction into the Kondo effect in the lattice is 
an outstanding problem, on which the impact of the insight from the two-impurity Kondo 
model has to be fully realized. 
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APPENDIX A: DERIVATION OF (H) AND (|2|) 

Under the transformation UHU~ l , only two terms in H are affected. They are 
h u /(2n)JZ dxd x ^ s (x) and 

= y JZ dx { n '( x ) + [^(*)] 2 } • ( A1 ) 

As for Hi, the transformation affects the term containing d x <& s (0)S+ apart from eliminating 
cos$ s (0) and sin$ s (0). Using the mode expansions, 
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dp 



UJx) 



2n^2\p\ 
dp\p\ 



2W2|p| 



-# s (p)e^ + ^(p)e"^' 



a|p|/2 



-a|p|/2 



(A2) 

(A3) 



we can write 



# w =^ n ^ip^mmp), 

J-oo Z7T 



d x $ s (x) = % I dpJ^e~ ap/2 UpV PX ~ <Pl(p)e~ vpx 
Jo V 2tt l 



P_ -qp/2 

The commutation relation for the Fourier components is 



(A4) 
(A5) 



[0.(p),0j(p')]=27rtf( P - I /). 



Next, let's introduce a generalized transformation operator 



(A6) 



17(A) = e - iXS +^(°) = e" iA5 + /r dpe ~ ap/2 [^) +< ^A 



r 2wp 



(A7) 



and define two A-dependent functions, 



h{\) = U{\)Hi s) U-\\), 
/ 2 (A) = ?7(A)^$ s (x)f>- 1 (A). 



(A8) 
(A9) 



We note U(X = 1) = U. Using the commutation relation ( |A6| ) , it is straightforward to verify 

d 2 



dX 



:A(A) 



a 



d r°° 

— / a (A) = -2/ dpe- Qp cos(ps)S* 
aA jo 



(A10) 
(All) 
(A12) 



From (|A10|) and ( |A11| ), we obtain 



UH^U- 1 = H$ - v F d x $ s (0)S z + + -(S z + ) 2 . 



a. 



(A13) 



From (|A12[) , we obtain 

Ud x $ s {x)U- 1 = d x $ s (x) - 2 / dpe~ ap cos(px)Sl. 

Jo 
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(A14) 



This implies 

2 



Ud^iOp- 1 = d x $ s (0) - ^S z + , (A15) 

A /"CO A /"CO 

E7 / dxdJisip^U- 1 = / dxd x $ a (x) - 2nS z + . (A16) 
Substituting ( gig ), (|AT5|) and flSIEQ into UHIJ- 1 , we obtain the results (gj) and (g2j) . 



APPENDIX B: DERIVATION OF THE EFFECTIVE HAMILTONIAN IN THE 

SECOND ORDER PROJECTION 

In this appendix, we shall derive the coefficients go, gi, V, a u and a s in the effective 
Hamiltonian ([14]), from the second order projection (j42[) . The Hamiltonian to be projected 
is © + © + ©. 



The first order contributions to the effective Hamiltonian are 

QHQ = -iv F r dx ]T ^\{x)d x ^ x {x) + i^L U sf (0) - ^(0)1 QSIQ 

/oo 
dxil>l(x)i/) 8 (x) + H 2 . (Bl) 
-oo 

In terms of the local fermion operators d and d\ we have QS y _Q = i(d — d>). Performing the 
transformation fl43|) to install the anticommutation relations between the extended and local 
fermion operators, the hybridization term between the extended and local fermion operators 
in (]Bl| ) becomes, 

-^$ sf (0)-$ f (0)](d+J). (B2) 

In the second order of the projection (fE|), the local impurity spin state is virtually excited 
from one of the doublet to either (| ||> — | |j>)/v2 or (| f|> +| ||>)/v / 2, then returns 
back to the doublet. From Figure the mixing terms between the doublet and the excited 
states are those in fl36|) + ( p7| ) which contain and S±. If the local impurity spins leave 
and return to the same state of the doublet, the generated contributions to the effective 
Hamiltonian have the form of either d^d or dd) . These are simply the renormalizations to 
the RKKY interaction, shifting its critical value determined by 
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K Z = K ± + 



(B3) 



where the omitted extra terms stand for the above mentioned renormalizations. The relevant 
term in the effective Hamiltonian around the critical point is, like the mass term in the usual 
critical phenomenon, 

K z + K x 



K c ) Sd, with K c = — (Jl + 



(B4) 



The other contributions come from the situations when the local impurity spins start 
from one state but return to the other state of the doublet. From Figure |3|, we see that these 
contributions must come from the projection of the product of S\ and S*L, 



Q 



Vl/W./(o)-V./(o)vI/(o) 



Sz 



X- 



l-Q 



-(K ± + T K ) + iv F IZo dx^ sf (x)d x ^ sf (x) - V V,/(0) - Vl/(0) ^/(0) + ^J(O) 



x 



v F Jj 



V2 



ira 



iP sf (0)+4 f (0)]S* + \Q + Q 



v F Jj 



l-Q 



W(0)+^I/(0) 



SI 



-(K ± + T K ) + iv F fZc dx4 f (x)d x ^ sf (x) - V [^/(0) - Vj/(0)] [^/(0) + $(0) 

^1,(0)^/(0) - ^-/(0)Vl/(0)] S*_] Q. (B5) 

Besides the energy gap K± + Tk between the doublet and the local excited state (| tt> 
+ | l|>)/v / 2 ! we have also kept the intraband terms in the intermediate denominator which 
will be expanded as 

(K ± + T K ) + iv F f°° dx^ sf {x)d x ^ sf {x) - V k/(0) - 4 f (0)] U/(0) + ^(0) 



IV p 



+ 



K ± + T K (K± + Tk) 2 J- 
V 



dxipl f (x)d x ip sf (x) 



Substituting flB6|) into flB5|), we obtain three contributions to the effective Hamiltonian. 
The first contribution is 



(B6) 



2(K 
+ 



T, 



K) 



P== { [vv(o) + <(o)] [ii f W. f (Q) - Mo)iif{o)] QS x + s*_Q 
^1/(0)^/(0) - V./ W£,(o)] [1M0) + vI/(o)] QS^q} . (B7) 
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Using the fact QS^S+Q = <$ and QS+S Z _Q = d, and carrying out the transformation ([13 
we simplify (p7| ) to 



u 2 T 1 - T z 



2{K ± +T K ){2^af/ 2 



^/(O)-^ (0)1 (d + d*). 



(B8) 



In deriving ( |B8|) , we have used the relation 

{^ / (0),^ t / (0)} = 5(0) = -^ 



27TOJ 



(B9) 



Noting that ( |B8| ) is a renormalization to (|B2|) . 

The second contribution to the effective Hamiltonian from 



is 



„.3 j± jz 
F + J n 



^ s /(0)+ Vl/(0)1 f° dx^ af {x)d x j, aJ {x) 



2(K ± + T K ) 2 V2 

V s t / (0)V s /(0) - -0 s /(O)-0l/(O)l g5^5*Q + [^(0)^/(0) - WJ/(0) 



X 



x / dx^lfWd^.fix) ^/(0)+^/(0) QS Z _S*Q\. 



(BIO) 



Commuting all fermion operators evaluated at a; = to one side and simplifying the products 
using anticommutation relations, we find that (|B10|) contains a term 



7,3 T± JZ 



8(K ± + T K y(27ia) 3 / 2 



M0)-$U0)] (d-d*). 



(Bll) 



This is the leading irrelevant operator and has dimension 3/2. Noting that the combination 
of local fermion operators appearing in (pil| ) is d — d\ not d + gF! This is a vital difference. 
The third contribution to the effective Hamiltonian from (|B5|) is 



2(K ± +T K ) 2 V2 



na 



{ [WW + ^1/(0)] [^/(0) - ^1/(0)] [^/(0) + v}(o) 



X 



^1/(0)^/(0) - lM°)V>J/(0)] QSIS1Q + [^(0)^/(0) - ^(0)^/(0); 

> a/ (0) - <(0)] [^(0) + ^}(0)] [V s /(0) + ^1/(0)] QS^-<?} . (B12) 



This contribution can be reduced to 



Vd 2 T 1 - T z 
v V F J + J m 



2(K ± +T K ) 2 (27ra) 5 / 2 



^(0)+ 4(0)1 (d-d*). 



(B13) 



35 



Again, we note that it is d — d) appearing in ( P13|) ! This is the second relevant operator 
which is present only when the particle-hole symmetry is broken, i. e. when V ^ 0. 
The staggered magnetic field coupling term comes from 



+ Q 



v F Ji 



> \J1i\a 
This term is simplified to 



'V./(o)+VI/(o) 



>./(o) + vI/(o)]^|q 

1-Q 



-(K ± +T K 



■{h.St}Q. 



(B14) 



h s v F Jj 



\p sf (o)+4 f (o)} {d-S). 



(B15) 



(Kl +T x ) v / 2^ 

To obtain the uniform magnetic field coupling term, we restore the bose field 4> s (x) 
through ipl(x)ip s (x) = d x $ s (x)/(2w) in (g6|) and (|3?p. Then we integrate out <p s {x) exactly. 
This is carried out as follows. First, we notice that the terms in ([H]) and fl3"2] ) containing <p s 
can be rewritten as, upon inserting d x Q s {x) = y/n [d x (j) s {x) — IT s (a;)], 



^ [°° dx \u 2 s (x) + [d x Ux)f + —d x $ s (x) + l±6(x)d x $ a (x)S 



^ r dx 



2 J- 



UJx) 



h, 



+ 



-i 2 



d x 6 s (x) + 



+ ?±6(x)d x $ s {x)S 

7T 



+ I > 



(B16) 



up to an additive constant. By introducing 



s (x) = s (x) + 



2^pnvp 



x. 



(B17) 



we can recast (P16| ) in the form 

v 4- r dxln 2 (x) + \djjx] 



2 J- 



2 + ?±5(x)d x $ s (x)Sl - ^±S z Mx) 

7T TIVf 



(B18) 



where and d x <& s (x) are correspondingly defined as, in consistency with fl2ll) and the 

relation fl s = d t <f) s , 



UJx) = UJx) - 



h, 



d x $ s (x) = d x $ s (x) + ^. 

V F 



(B19) 
(B20) 
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Since the uniform field h u only appears in the last term of ( P18 ), we only need to project it 
onto the lowest doublet in the next step. The contribution is 

r v F .F 



Q 



iMO) + 4/(0)] [4(0) - 4(0)] S z + + v F J z + i;t(0)il> s (0)S* 

1 -Q 



x LJ11-!^s' + \q+q\-!^I±s. 



X 



-Tk 



\fj sf (0) + 4/(0)] U/(0) - 4(0)1 s z + + v F Jtyt(0)ip a (p)S K +) Q. 



With a little algebra, one can show that (B2 1|) contains 



(B21) 



h u VFJ+Jt 



8ttT, 



K 



^ s/ (o) + 4 / (o)l 1-0/(0) -4(o) 



+ ^f^4(o)^(o). 



(B22) 



Combining the results (|E2|), (|B|), flBTT] ), QBEj), flBTgp , (jggg ) together and omitting the 
tilde signs on ?/>'s, we obtain the effective Hamiltonian (P>|), ( p)|) and (0). The qualitative 
physics does not depend on the numerical values of the coefficients go, g%, V , a u and a s in the 
effective Hamiltonian. It should be kept in mind that the obtained expressions, ( f£8|) to (p^D, 
for the numerical coefficients of the effective Hamiltonian should not be taken too literally. 
The purpose of this appendix is to illustrate how each term in the effective Hamiltonian 
arises from the projection rather than accurately determining the coefficients of the effective 
Hamiltonian. A practical way to determine them probably is to fit numerical results or 
experimental data. 



APPENDIX C: EFFECT OF THE MARGINAL OPERATORS (|55|) 



In this appendix, we shall show that the only effect of including the marginal particle- 
hole symmetry breaking operators (|55f) is to slightly renormalize the two basics energy scales 
Tk and T c . 

The marginal operators flSop correspond to the following terms in the action, 



/OO 
^—a f (uj n ,k) 

Our task now is to diagonalize 



a v a(-uj n ) - 



dk', 
2tt 



-u n , k') 



(CI) 
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S(b sf , a f ) = S 2 {b sf ) + S 3 {a f ) + 



(C2) 



where «S 2 (6 S /) and S 3 (aj) are given by (|70"D and ( |7T]) respectively. More specifically, we need 
to find new linear transformations other than (|76|) and ( |77D such that the hybridizing terms 
in (|C2|) are canceled out. The desired transformations are 



b sf (uj n , k) = b sf (uj n , k) + £i(a;„, k) a(u n ) + £ 2 (v n , k) b(u n 
a f (u n , k) = a f (u n , k) + £ 3 (u; n , fc) a(w„) + £ 4 (^n, k) b(u n 

where the four transformation coefficients are given by 

2iv F g + sgnu n a v V 2 /v 2 F 



(C3) 
(C4) 



fi(w n ,fc) 

fc) 

fc) 



iu n — v F k l + V 2 /v F 
2v F 



iu n — vpk 



gik 9 ; TAO ( V 1 \UJn\ 



v 2 F + V 2 
2iV a v - g sgnu n 



7T 



V F 



iuj n — v F k l + V 2 /vp 

2i V- 2VgiA/-K + Vg x \u n \/v F 



(C5) 
(C6) 
(C7) 
(C8) 



iu n — Vpk 1 + V 2 /Vp 

The ultraviolet cutoff A enters the transformation coefficients through the integrals ( |TDj ) and 
(pOl). In terms of the shifted Grassmann variables, the action (|C2T) becomes 



/oo dfc r \ _ ^ „ „ 

— S r(iw n - vf/c) a f (-u n , -k)a f (u n , k) + b sf (-u n , -k)b sf (u n , k) 
-oo Z7T L ^ 



— 5/(a; n , k)b sf (-uj n , k') > + 5 9en (a, 6). 



(C9) 



The generated local terms are 



[A^aa(|^n|)a(-^n)a(^n) + M b b( \u n \ )b(- UJ n )b(uJ n )} 



+ M ab (cf n ) a ( -a> n ) 6 (u n ) } 



(CIO) 



with 



2w F 



1 + V 2 /v 2 F 
2 

1 + V 2 /v F \v F 



21 



(Cll) 



7T / 
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2g 2 A 



V 2 \ 2 gi VV 



7T 



91 A 



(C12) 



1 + V 2 /v 2 - 



VV 2v F g l A 
9o 1 SlWn 



+sgnu; r 



F 

a v V 



71 



~ 2V 9l A V 9l 
V 1 



(C13) 



v F \ vr 'Up 

To obtain the total effective local action, we add the generated terms QUI 0| ) to (|73|). Intro- 
ducing a new local Grassmann variable, 



b(u n ) = b{u n ) + 



a v V 



-a(u n ), 



(C14) 



V - 2V gi A/7r 

we can compactly write the effective local action as, upon neglecting u> 2 terms in the matrix 
elements, 



Sill = -«E (a(-w„),6(-w, 

n>0 



X 



/ uj n /Z a + 2v F gl 



V 



5if + 2g g 1 uj n - r]uj r< 
The renormalized parameters are, 



-\5K + 2g g 1 uj n ) - r]u n 
u n /Z b + 2V' 2 /v F J 



a(u n ) 



(C15) 



go = 

9i = 
V> = 



go 



1 + V 2 /v 2 F 

gi 

V - 2V gi A/it 



5K = 5K- 



1 + V 2 /v 2 F ' 

2v F g (VV 2 9l A" 



1 + V 2 /V 2 F \ V 2 F 7T 



— = 1 + 



Tr 1+ 



y-2Vg 1 A/n / 



1 + 



Ag 2 A" 



V 2 /v 2 



1 T V - 1 Up 

a v V 



2g 2 A ( V 2 \ , 20x 



7T 



2(2^A/7r + ^VV/^) 
' V 2 /^ " 



(C16) 
(C17) 

(C18) 

(C19) 
(C20) 
(C21) 
(C22) 



V V - 2Vg 1 A/n L , . , . , 

effective action (|C15|) has essentially the same form as fl81| ) except for a wave function 
renormalization factor Z a and a new type term, f]u n , in the off-diagonal matrix elements. 



The 
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However, this new type term is irrelevant since it can only generate a (f]ui n ) 2 term in physical 
quantities such as free energy. At this point, it is clear that including the marginal particle- 
hole symmetry breaking terms with the coefficients V and a v V will only renormalize the 
two energy scales, the Kondo temperature Tk and the crossover temperature T c . 



APPENDIX D: CONTRIBUTION OF THE MARGINAL OPERATORS (g) TO 

THE STAGGERED SUSCEPTIBILITY 



In this appendix, we show that the contributions to the staggered susceptibility from the 



marginal operators (54) are negligible. 

The marginal operators ( |54"D have the following corresponding terms in the action 



Si 



stag 



roo dk 

2%h s 1^ / TT^fi-^n, k) 

„ J-oo Z7T 



/oo A^j 
—b sf (u nj k') + a"a(uj n ) 
-oo ZlT 



(Dl) 



Combining the last expression with the staggered field coupling term in (|69| ) and inserting 



the transformation (]76), we can write the complete staggered field coupling terms in the 
following form 



'stag 



2ih s J2 I ^a sf (-u n ,k) 

„ J — oo ZlT 



G' 



befiuJr,, k') 

oo 27T A ' ' 



+ (a s + a s g sgnu; n ) a{uj n ) + [a s H b{u n ) 



(D2) 



7T Vf 

The staggered susceptibility is obtained by calculating the second order perturbation of 



S 



stagy 



d 2 



dhi 



i 



q (S S t a gS S f a g} 



(D3) 



where the average is weighted by an action consisting of the free and decoupled Grassmann 
variables a s f, b s f and the effective local action (|54). The first term in S stag is a potential 
scattering term and does not mix with the other terms of S sta g in the second order pertur- 
bation. It thus gives a finite contribution to the staggered susceptibility and can be treated 
separately. Carrying out the calculation for (|T33|), we find that the singular part of the 
staggered susceptibility is 
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Xs = — 

vf 



2 



a s ^ -tfiA -^isgnw n (b{-u n )b(uj n )). 




n 



(D4) 



The propagator is given by, from (|84|). 



(b(-u n )b(u n )) 



Z b (\uj n \ + T K ) 



(D5) 



= —i sgnw. 



n 



•I + T c 2 + \uu\Tk 



where we have taken Ivpoi — Tr m the numerator of ( |D5| ) for simplicity. The singular- 



ity of the staggered susceptibility comes from the fact that at the critical point, T c = 0, 
(b(—u n )b(u n )) ~ l/{iuj n ) which gives rise to the logarithmic singularity for the Matsubara 
frequency summation in (|D4|) . From (P4j) , we see that the only effect of the marginal op- 
erators (|54D is to shift a s to a s + 2a' s giA/ir. Actually, one should be able to see this from 
( |D2| ) without doing calculation. 
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TABLES 





TABLE I. Definition of frequently used parameters 




symbol 


Definition (Eq. No.) 


Symbol 


Definition (Eq. JNo.j 


Dp 


Fermi velocity 


51 




PF 


Density of states 


V 




JZ T_L 


Kondo coupling constants 


a u 


if°il y oi u 




RKKY interaction 






Q x,y,z 


© 


5K 


(E3J 


h u , h s 


(H) 


A 




7-2,-L JZ,± 
J ± i J m 


(0) 


z b 


dD 


V 


P) 


T K 


dH) 






T c 


© 


90 
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TABLE II. The building blocks for constructing operators around the critical point. 



Operator 


Dimension 


Parity 


Particle-hole 


7r x-axis 


W(o)-<(o) 


1/2 


— 


+ 


- 


if>sf{0)+4 f {0) 


1/2 


+ 


+ 


+ 


V/(0)+^(0) 


1/2 


— 


— 


- 


^/(O) - v}(o) 


1/2 


+ 


+ 




d + S 


1/2 




+ 




d-S 







+ 






1 


+ 


+ 




d x 


1 


+ 


+ 


+ 
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TABLE III. All dimension 1/2 operators. The first one is the relevant operator. The second 
one could couple to the staggered field. The third one could become the second relevant operator 
if the particle-hole symmetry is broken. As explained in the text, [ip s f(0) — ijj s ^(0)](d — S) does 
not exist. This is because [tp s f(0) — V^j(O)] could hybridize with either d + d) or d — $, but not 
both. 



Operator Parity Particle-hole tt x-axis 

+ + + 

- + - 
+ - + 

- + + 



(d + <fl)(d-cft) 

[ij sf (o) + 4 f (o)](d-S) 

h/7(0)+^(0)](d-dt) 
[^(0)-4(0)](d-dt) 
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TABLE IV. All dimension 1 operators. The first one is the hybridization term. The second 
one can couple to the uniform magnetic field h u . The third and fourth operators could couple 
to the staggered field h s . The fifth and ninth are the marginal particle- hole symmetry breaking 
operators. 



Operator Parity Particle-hole tt x-axis 

ty sf (0)-4 f (0)](d + S) + + + 

[ij sf (o)-4 f (o)][ipsf(o)+4 f m + 

[Mo) + 4 f {o)](d + S) + 

[^/(0)-Vl/(0)][^/(0) + ^(0)] + - + 

hMo)- ^(0)1^/(0) -^(o)] - + + 

[^/(0) + ^(0)1^/(0) + ^(0)] - 

[V/(o) + v}(o)][V/(o)-v}(o)] - - + 

[^(0)+^(0)](d + cZt) + - + 

[V/(o)-v}(o)](d + dt) - + + 
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TABLE V. All dimension 3/2 operators. The first one is the only allowed leading irrelevant 
operator in the presence of particle-hole symmetry. 



Operator Parity Particle-hole n x-axis 

d x [i> sf (0)-4 f (0)](d-J) + + + 

3 x [V>/(0)+4(0)](d-dt) + + 

x [^(O)-^(O)](d-dt) - + + 

[4(o)Mo)-Mo)4m[Mo)-^l f m(d-S) + + - 

[^(0)^(0)-^(0)^t(0)][^ / (0)+^ / ( )](d-dt) + + 

[^(0)^(0)-^(0)Vj(0)][V-/(0)+^(0)](d-dt) + 

[^(0)^(0)-^(0)Vj(0)][^(0)-V}(0)](d-^) + 

[^t(0)^(0)-^(0)^t(0)](d + d t )(d-d t ) + + - 

hMO) - vI / (o)][^ / (o) + vI/(o)][V/(o) + $(o)](d - dt) - 

hMO) - vI/(o)][Vs/(o) + <(o)][V/(o) - ^(0)](d - dt) + + 

hMO)- ^(0)1^/(0) + vl / (0)](d + dt)(d-dt) + 

hMO) + ^(0)1^/(0) + 4(0)][^(0) - ^(0)](d - dt) + 

hMo) + vI/(o)][^/(o) +v}(o)](d + dt)(d-dt) - 

[^/(0)+^(0)]^/(0)-^(0)](d + dt)(d-dt) + 

hMO) - ^(0)1^/(0) + ^(0)1^/(0) - ^(0)](d - dt) - - + 

hMO)- ^(0)1^/(0) +^t( )](d + dt)(d-dt) + - + 

hMO)- ^(0)1^(0) -^(o)](d + dt)(d-dt) + + 

[^/(0) + ^/(0)]^/(0) -4(0)](d + dt)(d-dt) + + 
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FIGURES 

FIG. 1. The phase diagram of the two- impurity Kondo model. V is the energy scale charac- 
terizing particle-hole symmetry breaking strength. K is the fully renormalized RKKY interaction. 
Tk is the Kondo temperature. Except at the critical point marked by the black dot, the low energy 
behavior is Fermi-liquid type everywhere. The shaded area is the region where our solution applies. 
The radius of the solution region is a fraction of Tk- 

FIG. 2. The energy level scheme of the four impurity spin states. The up and down spin 
states refer to the eigenstates of the operators Sf and Sf. At —K z = K±, the two levels 
(I tl> — I IT>)/v / 2 and (| jj> + | ||>)/\/2 become degenerate, forming a doublet. The super- 
ficial degeneracy between the doublet and (| f|> — | ||>)/v / 2 is lifted by the Kondo interaction 
term J~. 

FIG. 3. The four impurity-spin states and the impurity spin operators connecting them, 
= S£±S$ for X = x,y,z. 



FIG. 4. The crossover function for the specific heat, Eq. (|9q), for various values of T c /Tk- Tk 
is the Kondo temperature, and T c is the crossover temperature. 



FIG. 5. The crossover function for the entropy, Eq. (|92|), for various values of T c /Tk- 
FIG. 6. The crossover function for the staggered susceptibility, Eq. (97), for various values of 



T c /Tk, normalized to its value at Tk/2. 
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